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Abstract 

Recent experiments in a mixture of two hyperfine states of trapped Bose gases show behavior 
analogous to a spin-1/2 system, including transverse spin waves and other familiar Leggett-Rice- 
type effects. We have derived the kinetic equations applicable to these systems, including the spin 
dependence of interparticle interactions in the collision integral, and have solved for spin-wave 
frequencies and longitudinal and transverse diffusion constants in the Boltzmann limit. We find 
that, while the transverse and longitudinal collision times for trapped Fermi gases are identical, 
the Bose gas shows unusual diffusion anisotropy in both dipole and quadrupole modes. Moreover, 
the lack of spin isotropy in the interactions leads to the non-conservation of transverse spin, which 
in turn has novel effects on the hydrodynamic modes. 

PACS numbers: 03.75.Mn,05.30Jp,05.60.Gg,51.10.-Fy,67.20.-hk. 
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I. I. INTRODUCTION 

In JILA experiments, a mixture of two hyperfine states was found to segregate by 
species. The theoretical explanation for this behavior is based on the two states playing 
the role of a pseudo-spin- 1/2 system, having transverse spin waves. The theory of these new 
effects is based on old ideas of the transport properties of polarized homogeneous quantum 
gases of real spins, such as ^He gas and solutions of ^He in liquid '^He, j9|,[iy] transcribed to 
the trapped gas pseudo-spin case. 

Besides spin waves, the theory for homogeneous polarized fermions or bosons led to 
the prediction of anisotropic spin diffusion in the degenerate state.P],[l^,|lll When a spin 
nonuniformity is longitudinal, that is, with a variation in the magnitude of the magneti- 
zation, the spin diffusion coefficient is D\\. On the other hand, in a spin-echo experiment, 
the magnitude of the magnetization is uniform but it varies spatially in direction. The 
corresponding diffusion coefficient, is less than D\\ when the system is polarized and 
degenerate. Experimentally this feature has been seen, but was not always in reasonable 



accord with theory. 



Further, Fomin 



12j has suggested the effect should not exist. However, 



a recent experiment j9j] has overcome several possible experimental objections and finds good 
agreement with theory. Moreover recent papers have presented theoretical analyses 
that question the validity of Fomin's approach. 

Thus it seems useful to see whether a similar difference between longitudinal and trans- 
verse diffusion in trapped gases might provide an alternative testing ground for this question. 
However, what we show here is that the physical possibility of having differing interac- 
tion parameters between up-up, down-down, and up-down states (interaction anisotropy) 
provides a new physical basis for anisotropic spin diffusion for bosons even in the Boltz- 
mann limit. [l^ For longitudinal diffusion in the Boltzmann limit only up-down scattering 
contributes. However, in the transverse case, two spins at differing angles approach one 
another, and the scattering can be analyzed as being a superposition of, say, up-up and up- 
down scattering. In the fermion s-wave case, the up-up part gives no contribution, and, in 
the Boltzmann limit, the diffusion coefficients are identical. In that case one must go to the 
degenerate limit to see the anisotropy, which then is expected to arise because the density 
of scattering states differs in longitudinal and transverse cases. |Q] On the other hand, for 
bosons, for which both the up-up and down-down scattering rates do contribute, we find 
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an anisotropy even in the Boltzmann limit, but only if the various scattering lengths differ. 
We have here the interesting effect that, although both gases obey Boltzmann statistics, 
there is a macroscopic difference between fermion and boson behavior because of quantum 
statistical effects in the two-particle scattering. Moreover in the Bose case a new element 
enters: the interaction now introduces a spin-spin relaxation mechanism (a T2 process) so 
that transverse magnetization is no longer conserved, j?!] For very small relaxation times r, 
the hydrodynamic relaxation rate is no longer is linear in r ; now it diverges. For fermions 
it remains well behaved. 

Below we first use the moments method to compute the spectra of the lowest-lying lon- 
gitudinal and transverse modes. However, with that method we obtain a transverse decay 
rate 7_l that diverges as r approaches zero, in contrast to the usual diffusive behavior where 
7_L oc r. In this case it is necessary to solve the local hydrodynamic equations to find the cor- 
rect behavior, in which the hydrodynamic solutions are localized at the low-density regions 
at the edges of the cloud where the collision time is longer. The result is a much smaller 
decay rate than that obtained with the moments method. 

In helium the dipole modes were detected experimentally, but in the trapped gases it was 
more convenient to look at the quadrupole modes. Q], ^Thus we calculate both here. We 
have previously reported results on the dipolar modes. 15| 

In our previous work, Ref . we derived an analog of the Landau-Silin equation for a 2 x 2 
density operator rip (here acting in the pseudo-spin space), with effective mean- field single 
particle energy matrix ip. We can write the density and single-particle energy in a Pauli 
representation as rip = ^ (^fpl + nip ■ and Sp = {^Cpl + hp ■ where 5" is a Pauli matrix, 
\{fp ± rnpz) give the diagonal components of the density ripi = Upn, and nip represents the 
polarization, which in equilibrium is along the axis z. We find the following approximate 
equation for nip: 

Pi dnip dU dnip 



drtip 2 ^ 
^--hxmp + E 



with mpz{r) = ripi — np2 and npi2(r) = n*2i = |mp_(r) = ^{rrip^ — irripy). The 2x2 collision 
integral is Ip. The effective mean magnetic field is 

h = ^z + r]—M (2) 

where 

= Vi-V2+ [(til - tu) rii - {t22 - tu) ns] (1 + r/) ■ (3) 



In these 77 is 1 (—1) for bosons (fermfons); M(r) = / dp/h^ mp(r); nj(r) = / dp/h^ np{r)] 
Vi is the external field for species i; U = |(Vi + V2); and = rii — n2. The fs can be 
evaluated in terms of the measured scattering lengths aap by using t^/? = 'ir^Tiaaii/m. 

The equilibrium solution in the Boltzmann limit is m^°) = exp[— /3(p^/2m + [/)] 

where N is the total number of particles, Ni is the number of species i, Ai = Ni — N2 is the 
total magnetization, and cu = (tUxCOyCUzY^^ . 

We have derived the collision integral for the Boltzmann case when the various interaction 
pa..a„,ente. diffe. agrees with the san,e quantity derived in Ke,.^.n,^ 

and reduces properly to previous results if all the t's are taken equal, [l^ , [l^ We find 



TT 



{a\Ip\a') = ^ I '^Pic?P2'^P3'^(Pi + P2 - P3 - P4)^(epi + ep2i 



-P31 



P2 /'a"0"2 l^Pl Ja-2(T 



~^a'a2 [('^Pi)o-ct'('^P2)o-20-2 + ''?('^Pi )o-CT2 ('^P2 )o-2cr'] 



"l"2t(7(72to-'(T2 [('^pa )o-0-' ('^p4)(T2 172 ~\~ fliP' 



P3 )o-a2 V^PA ) o-2cr' 



'1} 



(4) 



We will linearize the kinetic equation for nip around the global equilibrium value m^^^z 
and use a moment approach to compute the spin wave and diffusive damping just as done 
previously. The linearized longitudinal and transverse equations are 



dSm. 



pz 



dt 



E 



Pi dSrripz dU d5m 



pz 



m dvi dvi dpi 



^(T(a|Lp|a) 



(5) 



and 



d6m. 



+ t- 



+ E 



dt ' ^ h 
Pi d6mp+ 

m dvi 



dU ddrup 
dvi dpi 



2(2|Lp|l). 



(6) 



where Lp is the linearized form of Ip. 



In the next section we compute results for the monopole and dipole modes. Experiments 
have detected the quadrupole modes which we study in the third section. 
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II. II. DIPOLE MODES 



A. a. Longitudinal dipole. 

We use a variational function of the form 

Strips = (ao + aiz+a2pz)mf'^ (7) 

with X being position in units of {Pm)^^^'^ij^^ and pi momentum in units of {m/PY^'^. We 
take the 1, z, and Pz moments of the kinetic equation in both the longitudinal and transverse 
cases. The results for the longitudinal case, if we assume a time dependence of exp^iut) for 
ai and 02, are 

dao/dt = (8) 
iuai — ujza2 = (9) 
iuja2 + ujzai = -7»a2 (10) 



with 



and 



7„ = ft. (12) 

coming from integrating the collision integral. Eq. (jH)) indicates that the monopole mode does 
not decay in the longitudinal case, which is consistent with the conservation of magnetization. 
The effective external field Qo does not enter into the longitudinal modes. The second line 
is the magnetization equation of continuity. The relaxation rate 7|| agrees with that derived 
in Ref. 0. To find the modes one simply solves a quadratic equation. The dipole spectrum 
is plotted in Fig. ^as a function of t^^ = 1/7||, the spatially averaged collision time. In the 
small uJyT„ limit, one finds 



(13) 



which has the form of the lowest-order solution of a diffusion equation in a harmonic poten- 
tial. The other solution for small uJzT^^ is the high frequency mode 

uj = i—. (14) 
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The first of tliese solutions, tlie diffusive mode, is mainly a z mode while the second mainly 
a momentum dipole. For large oJzT^^ the solutions have roughly equal magnitudes of a\ and 
a2 with frequencies 

10 = ±uJz H • (15) 




0.0 0.5 1.0 1.5 2.0 

FIG. 1: Real (dashed) and imaginary (solid) components of longitudinal dipole spin wave modes 
versus average relaxation time ry. Note the linear dependence of Im(u;) for small ry characteristic 
of diffusive behavior. The dash-dotted lines represent the results of a numerical calculation to be 
discussed below. 



B. b. Transverse dipole. 

We again use the form of Eq. (0). We also assume that the effective external field Qq is 
quadratic in position, that is, the part Vi — V2 dominates the density dependent part in Eq. 
dni). We write 

no = A^(52 + f) + A,z2. 
Taking 1, z, and Pz moments of Eq. © yields the results 

dao/dt — iujQao = — 7rao (16) 
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i [u - uJo - i^W) ai - uJza2 = --7rai (17) 
i{uj — cum — uJo)a2 + uj^ai = —7^02 (18) 

where c^o and uj\\ are parts of the effective field: 



7T = 7o(l + r^)EPV^ 



trrrr — ^12 



!'12 
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with Z,^ = N„/N, and 



7± = 7„ 



7R - 3S 



(19) 

^I2 

withi?= (1 + r/) Ea^L/fT + (1-^7)^12 and 5 = 2ti2[(l+r7) ^(T<7/f7-'7^i2], and the mean-field 
frequency is 



where A is the thermal wavelength. Since ujq occurs in each equation we can remove it by 
including a time dependence e*'^°* in each af, that is, we go into the frame rotating at ujo- 
We assume this is done and all frequencies and will neglect this constant part. 
Comments: 

1) If the interactions parameters tij are all equal, we have 77^ = 0, i? = 5 = 2t'^ so that 
7_i_ = 7||. Eqs. (fT^ - (fTHj) then reduce to those of Ref. and the longitudinal and transverse 
relaxation rates are the same, which agrees with the standard result for a homogeneous real 
spin system in the Boltzmann limit. 

2) For fermions, we have 1] = —1, so that, even if the fs are not equal, 7t = and 
7± = 7||- In the s— wave limit only up-down scattering occurs for fermions and there can be 
no anisotropy in the Boltzmann limit. 

3) For bosons with unequal t's, the spatial averaged transverse relaxation rate is not 
generally the same as the longitudinal. Moreover, we have a T2-type relaxation rate for 
ao and in the equation of continuity (fTTj) . The interaction anisotropy behaves something 
like a dipole-dipole interaction allowing relaxation of the transverse spin, an effect noted 
previously in Ref. 3- 
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If, for now, we take 7t = 0, then the lowest mode in the hydrodynamic hmit {uJzT^^ ^ 1) 
takes the form 

0. = u;H + - (21) 

[l + {ujMr_L) J 

where t± = 1/7_l, and um = ujm — uj^^. The so-called "spin-rotation parameter" from Leggett- 
Rice systems is /i = uimTi/M- The form of Eq. (PT|) is the hydrodynamic frequency as 
modified by spin rotation. The first term in the fraction is the effective diffusion 

frequency while the second is the dipole-mode pseudo-spin-wave frequency. This diffusive 
mode is now shifted by the effective external field, uj\\. This mode is the z dipole mode; 
whereas the p-mode is found to be 

uj = um + — (22) 

which again diverges as t± gets small. 

For large uJzT± we find the pair of frequencies 

LJ = ±ujz + ^ {ujm + t^ii) + (23) 

The effect of non-zero 7^ is to allow a T2 relaxation of the transverse spins. The results 
are shown in Fig. 2, where we have set 1/t±_ = 7||, 77- = 0.02/r_L, and uTm = ^z- In the small 
T± limit, one no longer has the hydrodynamic decay rate approaching zero, but instead it 
diverges at the origin because Im(ci}) ~ (t^z^x +7t), and 7t ~ 1/t"±. However, although 
suitable for finite cur, the moments method is inadequate in the hydrodynamic limit. It 
fails because the simple forms assumed for spatial dependence cannot adjust to the spatially 
dependent relaxation rates. One must solve local equations numerically for the spatial 
behavior. 

To obtain the hydrodynamic equations we expand the momentum distribution in terms 
of Hermite polynomials 

5m,z = e-^P'/2„^ ^ ^^(^^ ^)^^(p) (24) 

fc=0 

Substituting this into the kinetic equations, neglecting terms in a;||, integrating over the 
momentum, and keeping terms lowest order in t±_ gives, in the transverse case. 



dt5M+ + dJ+ = --fT{z)6M+ (25) 
kT 

dtJ+ + — dJM+ + ujIz6M+ + iujM{z)J+ ^ ~l±{z)J+ (26) 



0.0 0.5 1.0 1.5 2.0 



FIG. 2: Real (dashed) and imaginary (solid) components of transverse dipole spin wave modes 
versus average relaxation time r_|_ when the transverse spin is not conserved. The dotted line 
shows the lower imaginary mode when the transverse decay rate 7t = 0. In this plot lom/^z = Ij 
u)\\ = 0. The linear behavior of lm.{uj) for small t_|_ characteristic of diffusive behavior is destroyed 
and replaced by a divergence within the moments method used here. 

where 5Mj^{z,t) = co{z,t) is the nonequilibrium magnetization density, J{z,t) = 
J dp/h^{p/m)6mpz = Ci(z, t) is the spin current, and7_L(2;) = 'yj_{0) exp^—jSmuj^z^ /2). Anal- 
ogous equations hold in the longitudinal case. On the RHS of Eq. (j^Uj) the k = 1 momentum 
distribution has been treated as an eigenfunction of the linearized collision integral. This is 
justified by a numerical calculation of the matrix elements of the collision integral, which 
gives 



The eigenvalues of the hydrodynamic equations have been calculated numerically for the 
dipole mode with boundary conditions 6M{0) = 0, J(0) = 1, and J{oo) = 0, and the 
monopole mode with boundary conditions 6M{0) = 1, J(0) = 0, and J(oo) = 0. For the 
longitudinal and isotropic transverse cases this leads to only small corrections to the r ^ 
part of the spectra obtained by the moments method as shown in Fig. d 



7±(0)(1.000i7i(p) + 0.123if3b) 




(27) 
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FIG. 3: Imaginary part of the spin-wave spectrum vs. 0JzT± for the monopole and dipole modes 
with ujm = ^zi It = 0.027^, oj\\ = 0, and 7^ 7|| for both the moments method (dashed) and 
hydrodynamic calculations (solid). 

However, for 77- > the hydrodynamic spectrum differs qualitatively from that of the 
moments method calculation. As the hydrodynamic dipole and monopole modes do 

not decay at a rate ~ 1/t_l, but instead decay at a slower rate. For tu^r^ <^ a = 'yT'T±, the 
hydrodynamic equations can be Taylor expanded about the center of an effective complex 



potential, giving a decay rate Uz^2a\og{y/a/LUzT±) (See Fig. IHl) At small enough uJzT± the 
T2 decay of the magnetization at the center of the trap causes the monopole and dipole 
modes to coalesce into spin-waves localized on the lower density regions on the left and right 
sides of the trap. (See Fig. HI) 

III. III. QUADRUPOLE MODES 

Experiments have investigated only the quadrupole modes. The increase in theoretical 
complication is considerable, especially when we include unequal scattering lengths. For the 
moments calculation we consider a function of the form 





for both longitudinal and transverse cases. 
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FIG. 4: Profiles |(5M(z)| of the the dipole modes vs. y/ (3mu-^z, normahzed with a Hermite weighting 
function. Tallest peak is the hydrodynamic mode for 77^ = and t± = 0.01. Middle peak is the 
moments method ansatz. Smallest peak is the hydrodynamic mode for jt = 0.027_|_ and cj^t = 0.01. 



A. a. Longitudinal quadrupole. 



The equations we find are as follows: 



ao 

iujai — uJzCi'2 
2..a, + ( + I) a. - 2..a3 
LJza2 + {iuj + 73)03 



-(ai + as) 









(28) 
(29) 
(30) 
(31) 



where 



73 = 711 



3 a '-12 



When all the t's are equal we have 73 97||/5 for bosons and these equations reduce to 
those of Ref. fl. 

We can again consider the solutions to the equations in limiting situations. For very large 
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cj^Tii (with Til = 1/711 ) we find the the three solutions 

(32) 



±2cj, + i(73 + 7|| 



,•73 
^ 2 



In the small ^zT^^ limit we again find a diffusive mode (proportional to a second order 
Hermite function in 2;^) : 

as well as two higher frequency decaying modes (one in zp^ and the other the Hermite in 

'"II . (33) 

^73 



Fig. gives shows the modes for all cj^r,. 



B. b. Transverse quadrupole 



Taking moments of the quadrupole magnetic equation Eq. (jH)) and again going into the 
frame rotating at ojq gives the following equations 



1 7 

iu (ao + ai + 03) — iu\\ai + 7t ( ao + + g^3 







It f 3 7 \ 
itu (ao + 3ai + 03) - iuon (oq + 602) - 2uza2 + — U^o + 7:01 + = 

2 V 2 o / 

iu (ao + ai + 803) - itUjiai + 26^^02 - 2icjAfO,3 + ^7t (^«o + ^^i^ + (303 = 



where 



C2 
C3 



370 1 + 



12 



12 



+ 6t 



^12 



- ^2 



12 



7o 



- + (1 + n) — 
3 ^ '^60 



/79 _25 



(34) 
(35) 
(36) 
(37) 

(38) 
(39) 



^12 a 6 

In the case of equal values of all the t's we have C2 7|| and C3 ^ 187||/5 for bosons. As in 
the dipole case the first equation represents a nonconserved monopole mode. If we take uj\\ 
and 7r to zero, this equation represents the zero frequency monopole mode. We see then 
that a non-zero uj\\ couples the monopole and quadrupole modes, because it arises from a 
quadratic external field. We will now consider solutions in various special cases. 
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FIG. 5: Imaginary (top) and real (bottom) components of longitudinal quadrupole spin wave modes 
versus average relaxation time r||. The diffusive mode is shown with darker lines. 

Zero values of 7^, ojm, o,nd uj\\: When all three of these parameters are set to zero, we 
get modes that look just exactly like the longitudinal case of Fig. El The diffusive mode's 
imaginary part goes to zero with uJzT2 as expected; here T2 = 1/C2- Two of the modes have 
degenerate relaxation times. 

Non-zero um] zero 77- and With non-zero mean-field value um the degeneracy is now 
split as shown in Fig.lHl For small uJzT2, the real parts of two of the modes have finite values, 
ujm and ujm/2,, but the imaginary parts of those modes still diverge so that they are very 
short-lived. The diffusive mode, whose imaginary and real parts both go to zero as T2 0, 
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is given, to order (c<j^r2)^, by 



00 



Au,T2{i - UJMT2) 



(40) 



This behavior is of the same form as the transverse dipole mode, Eq. ()21|) with the factor 
of 4 occurring for this higher order mode. The relaxation time T2 reduces to ry only when 
all the t's are equal or for fermions. So even ignoring the effect of 7t we would have boson 
diffusive anisotropy since T2 is not the same as t\\. The monopole mode decay rate vanishes 
in this case. 





FIG. 6: Imaginary (top) and real (bottom) components of longitudinal quadrupole spin wave modes 
versus average relaxation time T2 with ujm = 2wz, jt = 0, a;|| = 0. The degeneracy in some 
imaginary parts is lifted, and the diffusive mode (dark lines) now has a real part. 

Non-zero 7t; and zero uf As in the dipolar case the existence of the non-conservation 
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terms means a different behavior for small T2 as we now see in Fig. [7| From the moments 
calculation we find again that the imaginary parts of all modes, including both monopole 
and diffusive, now diverge proportional to 7^ ~ I/T2. The figure shows just the imaginary 
parts of the monopole and diffusive, indicating that the transverse magnetization is no longer 
conserved for bosons. However, as in the dipole case, a more accurate analysis shows that 



the small r behavior shows a divergence as y log (1/072X2). Fermions on the other hand 
behave as usual, because their scattering can depend only on ti2 and not on the difference 
in t values. 




FIG. 7: Imaginary components of transverse monopole and quadrupole boson spin wave modes 
versus average relaxation time T2 with 77^ = 0.002/t2 and = 0. The solid lines are the modes 
with non-zero and the dashed lines show the previous results with 77^ zero. Now these modes 
diverge at small r. 



Non-zero uj\\ ; zero ujm cind 77- : The quantity ul!\\ arises from the existence of an external 
quadratic potential difference between the two species. Because it is quadratic, it causes a 
coupling between the monopole and quadrupole modes. In Fig. |H1 we show the transverse 
quadrupole modes with a small value of > 0. We see that it does cause splitting between 
the previously degenerate relaxation times. Further, the real frequencies that vanished for 
small uJzT2 now are now finite and split with non-trivial structure. A particularly interesting 
feature for quite small uJzT2 is the near coalescing of the imaginary monopole and diffusive 
modes due to the quadratic coupling. The upper mixed mode mode, instead of going to zero 
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as 4:UJzT2, as it did when diffusive now has half that slope as demonstrated in Fig. 1^1 . 




0.0 0.2 0.4 0.6 0.8 1.0 



FIG. 8: Imaginary (top) and real (bottom) components of the transverse monopole and quadrupole 
boson spin wave modes versus average relaxation time T2 with = O.luJz- The darker lines are 
the diffusive modes and the dotted line is the monopole mode; these are now mixed at small T2- 
The very small T2 behavior of the diffusive relaxation frequency and all the real modes is altered 
by the presence of this coupling with the monopole mode. See the next figure. 

IV. IV. DISCUSSION 

We have shown that rather unusual differences between Bose and Fermi gases in the 
Boltzmann limit can arise in spin-diffusion experiments, when the scattering lengths for dif- 
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FIG. 9: Imaginary components of the transverse monopole (dotted line)and diffusive quadrupole 
(solid line) boson spin wave modes versus average relaxation time T2 with cj|| = O.lcj^. The dashed 
lines shows the uj\\ = case for comparison. 

ferent spin configurations are not the same, and when the effective field for the transverse 
modes does not vanish. For fermions there is no distinction between longitudinal and trans- 
verse modes. In the fermion case to see such an anisotropy one would have to go to the 
degenerate limit where phase space differences would lead to a difference between longitudi- 
nal and transverse modes. However, for bosons, longitudinal and transverse modes show 
striking differences when the scattering lengths are not the same. Indeed in the transverse 
case the magnetization is no longer conserved as it is in the longitudinal case; instead of a 
linear oJzT dependence of the spin diffusive mode one finds a divergence. 

In experiments on Rb, the interaction anisotropy is very small. To test the novel effects 
predicted here one might use Na, which has a difference in interaction paraments; Numer- 
ically we estimate that for ^^Na 7_l can differ from 7|j by as much as 14% with 7t/7± ~ 0.04. 
Interaction differences might also be induced by using Feshbach resonance methods. 

We thank Dr. Jean-Noel Fuchs and Prof. David Hall for useful discussions. 
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